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Abstract 



We study the effect of next-to- leading order contributions on the phenomenon of sym- 
metry non-restoration at high temperature in an 0(N{) x 0(N 2 ) model. 
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It is a well known fact that broken symmetries are usually restored at sufficiently high temperatures 
(l], |2|, [|. Nevertheless, as already noticed by Weinberg there are models, like the O(N) x O(N) scalar 
theory, for which the high temperature phase is less symmetric than the low temperature one in a certain 
region of the parameter space. 

As noticed by Mohapatra and Senjanovic Q|, the effects of symmetry non restoration can have inter- 
esting phenomenological implications in connection to the CP problem. The same phenomenon has also 
been used to provide non inflationary solutions to the monopole problem, as argued by Langacker and 
Pi ||] and more recently by Dvali et al. ||. 

Weinberg's analysis of the problem Q is based on a lowest order perturbative calculation of the 
thermal masses, the region of parameters giving symmetry non restoration being associated to the values 
of the coupling constants giving imaginary values for some thermal masses. 

In this letter, we will compute the next-to-leading order contributions to the thermal masses and 
analyze how the region of symmetry non-restoration gets modified by their inclusion. The reason for 
doing this is twofold: from one side, as these models contain more than one coupling constants, it may 
very well happen that next-to-leading order effects dominate over the leading ones, in certain regions of 
the parameter space. On the other side, when trying to achieve symmetry-non restoration in realistic 
models, such as those considered in ref. Q, one may have to take values of the coupling constants large 
enough for next-to- leading order corrections to become significant. 

For simplicity, we will discuss here only the simplest model presenting this phenomenon, that is, when 
the O(Ni) x 0(iV 2 ) symmetry is global. In this case, the theory is described by the renormalized euclidean 
Lagrangian 

L R = ^d IJt <p 1 d IJ ,(f>i-\-~d IJ/ <f)2d IJ ,(f>2 + V(<f>i,<f>2)+L ct , (1) 
where <fii( 2 ) is a iVi (N2) components real vector, V(<pi, </> 2 ) is a 0(N\) x 0(N2) symmetric potential: 
V(0i, <h) = \™\4>\4>\ + \™\^2 + ^Xi^iMif + ^A 2 /i 2£ (0 2 2 ) 2 + iA/i 2e (0 1( /) 1 )(0 2 </) 2 ) (2) 
and L ct stands for the counterterms Lagrangian, calculated in the MS scheme (/x is the dimensional 
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regularization scale parameter and e = -^)- As we shall see later, the crucial fact for symmetry non- 
restoration is that the mixing coupling constant A can take negative values. In order for the potential to 
be bounded from below, the coupling constants have to satisfy the following relation 

AiA 2 > 9A 2 . (3) 

We will show how next-to leading order corrections to the thermal masses are calculated, by first 
looking at the simpler 0(N) case, which can be obtained from the above Lagrangian by simply setting 
4>2 =0. 

L R = + ^m?0i0i + ^Ai/i 2£ (0i0i) 2 + L ct . (4) 

As explained in || , perturbation theory based on this Lagrangian breaks down at high temperatures 
as powers of the temperature can compensate for powers of the coupling constants making radiative 
corrections large. The remedy to this problem consists in a redefinition of the mass and the introduction 
of a compensating counterterm 

Lr = ^id^i + \m 2 ^\ + ^A^C&fc) 2 + L ct - ^a 2 Mi , (5) 

where 

M 2 = m 2 + cr 2 (6) 

and a 2 is a temperature dependent counterterm (not to be confused with the mass renormalization 
counterterm already contained in L ct ) which should be determined self consistently in such a way to 
eliminate from the self-energy the terms which diverge quadratically with the temperature. To illustrate 
how the mechanism works, it is sufficient to consider the case m 2 > 0. To one loop, one has for the self 
energy at zero momentum, 



Here h is defined as: 



1 r°° x 2 

h(y) = — 5- / dx ; . (8) 
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For small values of y, h{y) has the asymptotic expansion 
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h ^ = 2i-^-m^y ^t^-2) • (9) 

7 being the Euler constant. Notice that in eq.(Q) we have subtracted the following UV divergent part 
which is cancelled by the mass counterterm in the MS scheme: 

^ = -A^^±l[i + l- 7] . (10) 

From eq. (f?]) , one derives the " gap equation" : 

M 2 m 2 s (2 + Ni) (M\ 

-7^ = 7^ + ^— z—Lhl — ) , (11) 



rp2 rp2 11 g \ T 

which, after making use of eq. ([)]), becomes 

M 2 m 2 (2 + N 1 )fl M\ M 2 , M , , . 

2T = ^ + A 1 L__L2^___ j +0 (A 1 ^ln-) . (12) 

Solving this equation to lowest order gives the well known thermal mass, 

M 2 ^m 2 + ^±Il — \ 1 T 2 . (13) 
3 24 v ; 



This choice of M would in principle restore the validity of perturbation theory for — 3> 1 in the sense 
that higher loops corrections will be suppressed by powers of A with respect to the "tree level" thermal 
mass eq.(13). Notice that when eq.(fL3|) is introduced in eq.(^), one gets at high temperatures (T ^> m) 



^( 2 ^)T*. (14, 



7T V 3 24, 

This term is the next to leading order correction to the self-energy, and although it is not immediately 
obvious, does not get modified by higher loops corrections, as they are at least of order Af In Ai (see ref. 
@] for details). 

Notice that the same contribution could have been obtained more directly by just solving the gap 
equation to next to leading order, that is, by keeping the term linear in 



giving, for -jg > 1 



i^ = a±^M!_3(£±^v\v . (i6) 

O OA ~ \ O OA I \ > 
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At this stage, one would be tempted to improve the result eq.(16), by incorporating also the logarithmic 



terms in the expansion of h. This would give us a correction to ^ of order \\ In Ai, but, as contributions 
of this order arise also from two- loops diagrams, it would be necessary to include them as well for 
consistency, something that we will not attempt in this paper. 

As we said above, we were assuming that m 2 was positive. In fact, the case in which m 2 is negative can 
be treated in a completely analogous way, because one is working in the regime of very high temperatures 
where one " self consistently " assumes that the symmetry is restored. By self consistency we mean that 
one is able to find solutions to the gap equations of the unbroken phase giving a positive result for M 2 . 
In this respect, notice that the next to leading correction, although coming with the opposite sign than 
the leading one (see eq. (filf) , will not change the sign of the self-energy unless Ai becomes large, in which 
case any attempt of a perturbative calculation becomes meaningless. Thus, in this case the subleading 
correction cannot alter the symmetry breaking pattern at high T. 

We now turn to the 0{N\) x 0(N 2 ) model. There will be a set of two coupled gap equations, 

x\ = ^ + X 1 ^^h(x 1 ) + XN 2 h(x 2 ) , (17) 

TfP" O I AT 

x\ = ji + \2—^h(x 2 ) + \N 1 h(xi) , (18) 

where 

a* = Mt/T (19) 
Using the asymptotic expansion eq.(j|), and keeping up to the linear terms, eqs.([l7|-|l8|) become, for ^ » 1 

x\ = ci-Ai— — xx- X— x 2 , (20) 

Z47T 07T 

2 , 2 + 7V 2 Nt 
x 2 = c 2 -\ 2 — x 2 -X— xx , (21) 

Z47T 07T 

where we have introduced the constants c, as 

Cl = Al ^2- + A 24 ' (22) 



c 2 =A 2 ^^ + A- . (23) 
The constants c 2 ; give the thermal masses Mf to lowest order: 

Mf = c t T 2 . (24) 

As noticed in ref.pj, it is possible to take A negative and still satisfying the constraint eq.(||), in such 
a way that one of these combinations of coupling constants, say C2, is negative (notice that due to the 
constraint eq.(^|) ci is necessarily positive), namely: 

A<0 , (25) 

AiA 2 > 9A 2 , (26) 

In this case, to lowest order, eqs.( fL7| - |l8| ) have no self-consistent real solutions and the corresponding region 
of parameter space was identified by Weinberg as the region of symmetry non-restoration (or rather of 
symmetry breaking at high T, if the symmetry was unbroken at T = 0). The question that now arises is 
whether the subleading corrections represented by the linear terms in eqs.( p0| - |2~l| ) can alter this picture. 
Notice that, as opposed to the simpler O(N) case, we now have more than one coupling constants and 
thus next to leading order effects could become dominant even when all the couplings are of the same 
order of magnitude. Now, eqs.(|20|-pT|) represent a pair of parabolae in the plane X\ and Xi which intersect 
in the upper right plane (then giving an acceptable solution to the gap equations and signalling symmetry 
restoration at high T) if and only if: 



8^< -Ar 2 + iVl 



2 + iVA 2 

Ai + ci 

48tt j 



XNi " 48tt 

The corresponding region of symmetry non-restoration is then described by the set of inequalities: 



(28) 



A<0 , (29) 
AiA 2 > 9A 2 , (30) 
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|A| > A 2 



2 + N 2 
3Ni 




2 + Ni 
48tt 



+ Ui 



2 + Ni 
48tt 



) 



2 



+ C1 



2 



(31) 



By comparing with eqs.( |25||j7| ), it is clear that the inclusion of next-to- leading order corrections reduces 
the region of the parameter space in which symmetry non-restoration takes place. We observe that even 
if there is always a region of symmetry non restoration, its size becomes smaller and smaller as N\ and 
N2 grow large. 

We can think of our model as representing the Higgs sector of a SU (5) grand unified model, like the 
one considered in ref. ||. In this case, taking Ni — 90, (j>\ represents the real components of the Higgs 
in the 45 representation of the group SU(5), while taking N2 = 24, (j>2 corresponds to its adjoint 24 
representation. 

A qualitative indication of the importance of the subleading terms can be obtained by plotting the 
regions ( p5[p7| ), (|29]-|3l|) at fixed values of one of the coupling constants and of Ni and N^. Figures (1) and 
(2) display these regions for Ai =4/5 and A2 = 1 respectively and N\ = 90, N2 = 24, which correspond 
to typycal values of the ST/ (5) parameters considered in ref.S. It is clear from both figures that the 
subleading terms have a non negligible effect, as they reduce (in the considered region of parametrs) by 
a factor of roughly two the size of the region of symmetry non-restoration. 

In conclusion, we have shown that the inclusion of next-to-leading order contributions can modify in 
a substantial way the symmetry breaking pattern of an 0(N±) x 0{N2) model at large T. Naturally, the 
phenomenologically interesting models are those in which part of this symmetry is gauged. Consequently, 
our analysis of the gap equations should be generalized to include this case. The results presented here 
lead us to expect that subleading corrections will play an important role also in grand-unified models. 
We expect to report on this issue in a forthcoming publication. 
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FIGURE CAPTIONS 



Figure 1. [Region of symmetry non-restoration at fixed Ai = 4/5 for Ni = 90 and N 2 = 24. The zeroth 
order region is the one enclosed by the solid curves, while the region obtained by including next-to-leading 
order corrections is the one enclosed by the upper solid line and the dashed line.] 

Figure 2. [Region of symmetry non- restoration at fixed A2 = 1 for Ni = 90 and N 2 = 24. The zeroth 
order region is the one enclosed by the solid curves, while the region obtained by including next-to-leading 
order corrections is the one enclosed by the upper solid line and the dashed line.] 



8 



